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Abstract. The dynamics of non-polar diatomic molecules interacting with a far-detuned 
narrow-band laser field, that only may drive rotational transitions, is studied. The rotation 
of the molecule is considered both classically and quantum mechanically, providing links to 
features known from the heavy symmetric top. In particular, quantum decoherence in the 
molecular rotation, being induced by spontaneous Raman processes, is addressed. It is shown 
how this decoherence modifies the rotational dynamics in phase space. 



PACS numbers: 33.80.-b, 03.65.Yz, 42.50.Ct 

1. Introduction 

Molecules interacting with electromagnetic fields foster a variety of dynamical phenomena. 
Specific wave packets of the internuclear vibration have been excited by shaped fs laser 
pulses [1]. Such vibrational quantum states have then been reconstructed by molecular 
emission tomography [2,3]. Furthermore, control of the internal molecular quantum state 
could be achieved [4-8], for eventually enhancing chemical reactions. 

Besides these applications of pulsed fields, also interactions with cw laser fields or static 
fields reveal interesting effects. For instance, the axis of paramagnetic molecules can be 
aligned in pendular states by applying magnetic fields [9, 10]. Furthermore, it has been shown 
how electric fields align polar molecules [11]. Non-polar molecules, on the other hand, seem 
to be rather isolated from the influence of electric fields. They are infrared inactive and thus 
their rotation is essentially undamped and isolated from the electromagnetic vacuum. Clearly 
this does not hold for the internuclear vibration and, in addition, the rotational statistics 
determines via ro-vibrational coupling its relaxation and decoherence properties [12, 13]. 

Nevertheless, a permanent dipole moment being absent, non-polar molecules still can be 
polarised by electrical fields and thus feel an effective alignment force when an anisotropic 
polarisability prevails. The quantum-mechanical eigenstates and energies have been discussed 
in Ref. [14], where it has been shown that the dependence of the eigenenergies on the 
field intensity can be used for providing a dipole trapping force. Recent experiments have 
implemented such a dipole trap for molecules [15]. 

In general the topic of producing [16-21], cooling [22,23], and trapping of ultra-cold 
molecules [15, 24] became a major focus of research in recent years. Given the ultra- 
low temperatures of a condensed atomic or possibly molecular gas, the molecular rotation, 
representing the lowest energy scale apart from centre-of-mass motion, then should play an 
important role. For sufficiently cold molecules one might expect decoherence first to appear 
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Figure 1. Diatomic molecule in an electric field pointing in z direction. The angles i3 and Lp 
specify the orientation of the molecular axis with respect to the laboratory frame. 



on the rotational energy scale before affecting the energetically higher degrees of freedom, 
such as the internuclear vibration. Moreover, the molecular rotation provides a toolbox for 
implementing various textbook examples of mechanics, such as spherical pendula and the 
motion of rigid bodies. 

In particular, we will focus here on homonuclear, diatomic molecules in a far-detuned 
linear-polarised laser field. This is a setup relevant also for trapping of photo-associated 
ultra-cold molecules in a dipole trap [15]. In Sec.|2la brief review is given on the classical 
description of the motion of the molecular axis. A link to a quantum description in terms of 
stimulated two-photon Raman processes is then presented in Sec. |3l Moreover, the theory 
is extended there by including spontaneous Raman processes in form of a quantum master 
equation. The effects due to these spontaneous processes on the molecular rotation are 
evaluated in Sec. 0] where we provide a phase-space picture by using a Wigner function for 
the molecular rotation. 



2. Classical mechanical analogue 

2.1. Classical dynamics 

Consider a diatomic molecule in an electric field of amplitude i?o pointing towards the z 
direction. The potential energy in such a field is given by 

V {^)^-dEn cos d-lEl/^acos^d, (1) 

where d is the permanent dipole moment of the molecule and Aa = a\\ — a±_ denotes 
the anisotropy of the polarisability with respect to fields parallel and perpendicular to the 
molecular axis. The direction of this axis with respect to the electric field is specified by the 
enclosing angle d, cf. Fig.^ 

The molecular axis specifies the z'-axis of a body-fixed coordinate system, whose 
orientation with respect to the laboratory frame is specified by Euler angles. From the 
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Lagrange equations one obtains the equations of motion for these as 

e(j?-(^2sint9cost9) + ^^^ = 0, (2) 

if sin^i? = const, (3) 

where 9 is the moment of inertia perpendicular to the molecular axis. Moreover, the 
conservation of energy E reads 

e(i?^ + ^2 sin^t?) = 2 [E -¥{■&)] . (4) 
Note that for our case the relevant Euler angles d and (p coincide with the spherical angles as 
depicted in Fig.^ 

The angular-momentum L can be expressed in the laboratory frame by the angles and 
their time derivatives as 

(I? cos (p ~ (p sin cos -d sin ip \ 
-d sin if + (p sin cos i? cos v • (5) 
ip sin^i} j 

From Eqs (|3} and (|5} it is seen that the component of the angular momentum in direction of 
the applied field is a constant of motion: = const. Thus the tip of the angular-momentum 
moves in a plane normal to the electric-field direction. Using the constant of motion Lz, 
Eq. can be rewritten as 

^=i?-C/W, (6) 

where the new potential reveals a centrifugal-type barrier 

I± 
sin^?9 

For later purposes we have introduced here the potential depths Ud = di?o, Va = AaE'Q /2 
and the z-part of the kinetic energy T^^^L^j (20). 

A molecule with permanent dipole moment, d^Q, omitting its anisotropic polarisability, 
an « a^, reveals a potential of a spherical pendulum, as shown in Fig.|2] The motion in ip 
with constant angular momentum is superimposed by a nutation in i9, which is equivalently 
described by a heavy symmetrical top with vanishing moment of inertia along the (molecular) 
body axis. The harmonically approximated nutation frequency is in this case (x)= \JUd/Q- 

However, for the case of non-polar (d — 0) but polarisable molecules (a|| > a±^, which 
is the focus of this paper, the anisotropic polarisability produces a cos^i^ potential, cf. Fig.|3] 
which for 7^ cannot be cast into the form of a spherical pendulum. Nevertheless, a 
spherical-pendulum-type motion can be expected, since also for this case the solution i?(t) can 
be obtained in terms of elliptic integrals. Expressed in terms of Jacobian elliptic functions [25] 
it reads. 



C/(i9) = - Ud cos -d-Ua cos^d. (7) 



cos i?o cn(ti;t, m) — (cos^i^q — Ci)(c„ — cos^zJq) Jti/cu sn{ujt, m) dn{wt, m) 
COSij[t) ~ - - 7 a \ 57 I ^ '/ ''^ 

1 — (Cm — cos^i/o) sn^(cjr, m) m/Cu 

where i?o is the initial angle at time t = Q. The constants appearing in Eq. (|8} are given as 

(9) 



e- 



-1' 



ci.u = -T\ -r^~i^z, (10) 

m V 4to^ 



2f/a 



Q m 

where the scaled energies e = E/Ua and Kz—T^/Ua have been used. 



(11) 
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Figure 2. Nutation potential U{'d) for a polar molecule scaled by the potential depth Uj^. SoUd 
(dashed) curve corresponds to Tz/Ud = (0.1). 
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Figure 3. Nutation potential C/(i?) for a nonpolar, polarisable molecule scaled by Ua- Solid 
(dashed) curve coiTesponds to ft^ =Tz/Ua = (0.1). 

2.2. Classical types of motion 

In the asymptotic limit of large energies e ^ 1, sn, cn sin, cos and dn ^ 1 so that 
cosi?(i) = cos{'dQ + ujt) results, which corresponds to the free rotator case with angular 
velocity uj ^J2E/Q. In this limit all three vector components of the angular momentum 
are conserved quantities. 

For lower energies the potential barriers lead to modifications, supporting several types 
of motion. For k, = unbound motion occurs for e > where the amplitude of the oscillation 
of d attains the full tt range, see solid curve in Fig.0] Given an initial angular momentum only 
in X direction, the molecular axis then continuously rotates around the x axis so that while L^. 
is modulated due to the potential barriers, it does not change its sign, see solid curve in Fig.|5] 
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Figure 4. Nutation angle i? over scaled time iot for Kz = 0, i9o = 0, ij;(0) = and 
IQlJal^ =0.025. Scaled energies and x-components of angular momentum are: e = 0.024, 
L2;(0)/n=0.16 (solid curve); e = -0.1, La;(0)/?i= 0.15 (dashed curve). 
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Figure 5. Lx/h over scaled time u)t for the same parameters as in Fig. 4. 

On the other hand, for the same case of = 0, for energies e < the angle d is subject 
to a reflection at the light-induced potential barrier. Thus the range of d values is restricted, as 
seen for the dashed line in Fig.^] In addition, since i} changes its sign at the reflecting barrier, 
the angular momentum correspondingly crosses zero and thus rotations with positive and 
negative orientations subsequently interchange, cf. dashed line in Fig.|5l Note that for > 
the motion of d is additionally bound by reflections at the centrifugal barriers near i9 = and 

TT. 

3. Quantum dynamics of the molecular rotation 

3.1. Implementation via two-photon Raman transitions 

For implementing the dynamics described so far, not only static electric fields can be used. 
Also linear polarised cw optical fields, that are far detuned from electronic resonances, have 
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Figure 6. Excitation scheme for implementing rotational Raman transitions. 



been shown to produce a similar interaction [14]. In this way one can take advantage of large 
polarisabilities of non-polar molecules to attain deep potential wells for the nutation angle. 
The specific excitation scheme for a dimer, depicting only the two lowest electronic potential 
surfaces that provide a dipole transition, is shown in Fig.|6] 

In order to avoid resonant vibronic transitions and to be left with the sought two-photon 
processes, a sufficiently large detuning A from the bare electronic resonance frequency is 
required. Given an optical field of spectral width much smaller than vibrational frequencies 
Wv, Stimulated Raman transitions occur that only affect the rotational degree of freedom. 
However, the relative rate of spontaneous vs stimulated Raman transitions, F/A, with F 
being the bare electronic linewidth of the dipole transition, shows that even for large detuning, 
A ^ F, a certain number of spontaneous processes occur These processes not only lead to an 
incoherent effect on the molecular rotation, such as relaxation and quantum decoherence, but 
also excite the internuclear vibration via Franck-Condon transitions. 

Nevertheless, as long as the laser remains detuned from resonant vibronic transitions, the 
electronic and vibrational degrees of freedom can be eliminated to obtain equations of motion 
for the molecular rotation alone. This is achieved by consistently eliminating the electronic 
coherences and electronic excited-state populations, keeping only terms up to second order 
in the molecule-light interaction. In Franck-Condon approximation the vibrational degree 
of freedom can then be traced over [26]. The requirement of avoiding resonant vibronic 
transitions leads to a maximum interaction time on which the molecule stays vibrationally 
sufficiently cold, which can be estimated as 



Here JIr = |(igo-E'oP/(4?i^A) is the Raman Rabi frequency with dg^ being the electric-dipole 
matrix element between electronic ground and excited states. The linear polarised electric 
field is specified by E{t) ~ Eq cos(tJoi) with ujq being an optical frequency far-detuned 
from vibronic resonance and fiz being the unit vector in z direction. Using these parameters 
the light-induced potential depth results as Ua = fi^R- 



(12) 
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The parameter 7y is the ratio of the difference of internuclear equiUbrium distances in 
ground and excited electronic states to the spatial extent of the vibrational ground-state wave 
packet. Since ry < 10 for alkali dimers, in the weak-field regime [$7^ < h/{2Q)] this time 
corresponds to a very large number of free rotational periods. 

Under these conditions the equations of motion for the orientation of the molecular axis 
in the far detuned optical field can be cast into a master equation of Lindblad form [26]: 



a = --[T + V,cr] 



J2 (s,aS}-^S}S,a-^aS}S, 



(13) 



i—x.y,z 

Here a is the reduced rotational density operator of the molecule and the free kinetic energy 
is given by [27] 

f = JV(2e). (14) 
The stimulated rotational Raman transitions then provide the sought potential in the form 



(15) 



where in terms of spherical angles the unit-vector operator reads if\n2 = cos-d {d, 
so that 



(16) 



recovers the potential as discussed for the classical treatment. 

The spontaneous rotational Raman transitions are given in terms of the operators Si, with 
i = x,y,z denoting the polarisation of the spontaneously emitted photon. They read 



S, 



Tor . . 
— :— rii n. 



(17) 



with fii (i = x,y, z) being unit vectors pointing into the corresponding directions. In the 
representation of spherical angles these jump operators can be cast in vector form to: 

cos ip sin cos 
sin (p sin cos 
cos^ 1) 



(18) 



3.2. Stimulated Raman processes 



For vanishing bare electronic linewidth, F = 0, a unitary evolution is recovered from the 
master equation il3\ that can be reformulated in terms of the Schrodinger equation 

ihdtl'f) ^if+V)\^), (19) 

with the Hamiltonian given by Eqs il4\ and ( I15> . In the \d, ip) representation the evolution of 
the wavefunction ^'(t?, (f) = {•&, 'p\'^) reads 



1 



sin i3 



(sin i9 dt) 



sin^^^' 



(20) 



Using the ansatz tp) ~ S'i,„(z) exp{im(p), with hm = Jz being a constant of motion, 
and z = cos the energy eigenstates are obtained from the angular oblate spheroidal wave 
equation [25] 



29 



9. [(i-^')a. 



+ Uaz'^ + Eir,, \ Slm{z) = 



(21) 
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Figure 7. Change of energy eigenvalues due to the Hght interaction for the value 2@Ua/h^ - 
1. 



The corresponding eigenenergies can be expanded in powers of the interaction potential 

Ua as 



26 



f4 
2 



1 + 



(22) 



(2Z-l)(2; + 3)_ 

They manifest a partial removal of the degeneracy with respect to m, as shown in Fig.0 The 
properties of these eigenstates and their angular confinement have been discussed in detail in 
Ref. [11]. 

The quantum dynamics of a system prepared quite analogously as discussed for the 
classical treatment, i.e. with an initial average angular momentum pointing only in x direction, 
is shown in Fig.|8] Here an initial coherent angular-momentum state [28] of the form 



|j, a, /3) 




[ sin(a/2)] [ cos(a /2)] e-^(3+-^)P 



(23) 



for j = 2, a = 7r/2, and (3 ~Q has been used. Classically that would correspond to the case 
K2=Q where reflections in the motion of -d can only occur due to the light-induced potential. 
But from Eq. (I23> it is clear that we deal with an initial state containing states \ j, m) with all 
possible values for m: 



1*1) = |j = 2,a = 7r/2,/? = 0) 

= i |2,-2) + i |2,-1) + yj |2,0) + i |2,1) + i |2,2). 



(24) 



Thus part of the wave packet in d will also be reflected by centrifugal barriers. 

In Fig.|8]an oscillation of ( J^) between positive and negative values is observed that is 
characteristic for reflections of the wavefunction at potential barriers. However, the transient 
behaviour at the crossing of zero indicates that part of the wave packet in is transmitted 
through the barrier, and that reflected and transmitted parts interfere for the short time where 
they spatially overlap. 
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Figure 8. Dynamics of {Jx{t)) over scaled, dimensionless time t = th/{20). Tlie laser 
interaction strength is Mlii/{2@) = 0.1 and the detuning F/A = 0.01 (solid line). The 
dotted line corresponds to F = 0. 

3.3. Spontaneous processes 

When spontaneous processes are included in the dynamics, i.e for F > 0, the full master 
equation il3\ must be solved. Due to the scaling of complexity of this numerical problem with 
the number of required rotational levels, it is advantageous to employ a quantum trajectory 
method. Then individual realisations of state vectors can be calculated and finally an ensemble 
average over the set of realisation reproduces the sought density matrix. 

Individual trajectories are comprised of non-unitary evolutions with the effective 
Hamiltonian 

HoS = f + v(^l + ^y (25) 

describing the conditioned evolution when no spontaneous processes occur, intermitted by 
spontaneous Raman processes described by the application of the jump operators Si. As can 
be seen from their definition ( I18> . these operators may excite the motion in the angle ip, by 
producing via Sx and Sy a weighting of the corresponding probability amplitude ^{"ff, <p) in 
(fi. Thus motion in Lp and consequently the angular momentum component are getting 
excited by spontaneous Raman processes. 

This can be observed in Fig.|9] where the variance of Jz is plotted over time. For F > a 
monotonous increase of the variance is observed, indicating that the kinetic energy is subject 
to heating. Of course the spontaneous Raman scattering will also lead to a suppression of 
coherent processes and the oscillations of the average angular momentum will be damped, cf. 
soUd curve in Fig.|8] 

4. Decohering dynamics in phase space 

4.1. Spherical Wigner functions 

In the same way as the angular-momentum coherent states are a generalisation of the 
familiar harmonic-oscillator coherent states, also the concept of phase-space quasi-probability 
distributions can be defined on the angular-momentum phase space [29]. Here we employ a 
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Figure 9. Variance of Jz over scaled, dimensionless time, t = th/ (2©). The pai'ameters are 
the same as in Fig.lsl 

Wigner function, which in the (2j + l)-dimensional Hilbert space of angular momentum j, is 
defined in terms of the density operator a and the spherical harmonics, Ysm(6', ip), as defined 
inRefs [30,31] 

WA^^ f) = \/^E E ^) ^^^Tlm^). (26) 

s— — s 

The multipole operators are defined by 

= V2iTT 5] (-1)^-"' ( ^ ^„ ) bm')(jm"|, (27) 

ni' 7n" 

with {Jll^^^^ being the Wigner 3j symbol. The factor in front of Eq. \2b\ ensures 
normalisation to unity. 

As the laser interaction will result in a change of the rotational quantum number j of the 
molecule we consider a sum of Wigner functions over all possible values of j: 

oo 

W^(i?,^) (28) 
i=o 

Albeit the sum over j, this phase-space distribution is not a complete description of 
the rotational quantum state. However, it proves useful for illustrating the laser-induced 
dynamics. For alternative approaches defining phase-space distributions independent of the 
particular value of j, see Refs [32-34]. Moreover, it should be noted that even though the 
angular-momentum coherent states are analogous to the harmonic-oscillator ones in many 
respects, the Wigner function of angular-momentum coherent states is weakly negative. This 
negativity decreases with increasing value for j and reaches zero only in the (classical) limit 

4.2. Time evolution of the Wigner function 

Figure [Tol shows the time evolution of the system starting initially from the state l^"!), 
cf. Eq. (I24> . The phase-space distribution is shown for a scaled time r = th/{2Q) = 660 
after the first full period of oscillation of (Jx) (upper left plot), and then at later times in steps 
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of quarter periods. One can see that the system almost returns to its initial coherent state after 
one period, when the spontaneous emission is neglected (left column). The peak in phase 
space has then returned to its initial position at ?9 = 7r/2 and (/? = 0. In the course of time the 
peak splits and after a half period, at r = 990, two peaks appear at i? = 7r/2, (^s = 0, it. That 
corresponds to an average angular momentum pointing now in negative x direction. Note that 
between half and full cycles, at r = 830 and r ~ 1165, the phase-space distribution develops 
several smaller peaks and also negative regions of the Wigner function appear. Including 
the spontaneous Raman processes, the structures in phase-space generally are smoothed and 
negativities are suppressed. This can be seen from the right column of plots in Fig. 1101 

As another example, in Fig. a superposition of two coherent angular- momentum 
states. 



has been taken as the initial quantum state. This state has a strong negative region of the 
Wigner function centred at = 7r/2, = 0. The plots of the Wigner function for this initial 
state show the same oscillatory behaviour as Fig.^| with the negative peak being shifted to 
= 0, TT at half periods. 

However, comparing the time evolution of these two cases one can see that the difference 
between the unitary and non-unitary evolution is much more pronounced in the case of the 
coherent state. Especially at the times between completed half and full cycles (at t = 830 
and r ~ 1165), when also the Wigner function of the time-evolved coherent state shows 
relative strong negative regions, one observes that these negativities have become weaker 
in the presence of spontaneous emission. The populations in the j-level manifolds for the 
final interaction time of Figs ^| and ^2 shown in Fig.^] It can be seen that for both 
the coherent state and the superposition of coherent states the populations in j = 0, 4 have 
increased by a few percent at the expense of the population in the initial manifold j — 2. 

That the coherent state in this example actually is much more sensitive to decoherence 
may seem counter-intuitive at first sight, but one should keep in mind that the angular- 
momentum coherent state is not classical in the usual sense of a harmonic-oscillator coherent 
state. That is, as has been seen above its Wigner function in general is not positive everywhere. 
Furthermore, angular-momentum coherent states apparently are not in general the most robust 
states against decoherence, as are harmonic oscillator coherent states for linearly coupled 
reservoirs. 

The substantial difference in how these two cases are affected by decoherence due to 
spontaneous processes can also be observed in Fig.^] where the time evolution of the purity 
of the density operator is shown. It can be seen that the superposition state M9\ develops 
much slower into a mixed state compared to the coherent state. 

The type of decoherence mechanism is in general determined by the form of the operators 
Si, see Eqs ( I17> and M%\ . In our case these operators differ from raising and lowering 
operators J± = Jx±iJy, and thus our decoherence mechanism itself is different from that 
discussed in Refs [35, 36], though for special choices of initial conditions particular features 
may be similar. 

5. Summary and conclusions 

We have considered a diatomic homonuclear molecule in a far-detuned, linear-polarised light 
field. First the classical motion of the molecule was studied. For the case of a molecule 





Figure 10. Time evolution of the colierent state | i ) . From top to bottom the figures show 
the Wigner function at scaled times t = 660, r = 830, r = 990, and r = 1165. Parameters are 
the same as in Fig. Island r = th/ {2&). 



with a permanent dipole, the potential was seen to be that of a spherical pendulum. For non- 
polar molecules, however, the potential is of a different form, but the solution for the nutation 
angle still could be found in terms of elliptic integrals. The motion was then treated quantum 
mechanically, where stimulated and spontaneous Raman processes were taken into account by 
a master equation. Solving this equation both for an angular-momentum coherent state and 
a superposition state, the effects of the decoherence induced by the spontaneous processes 
where compared. It was seen that, in this particular example, the superposition state was 
more robust to decoherence than the coherent state. 
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